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We investigate phantom cosmology in which the scale factor is a power law, and we use cos- 
mological observations from Cosmic Microwave Background (CMB), Baryon Acoustic Oscillations 
(BAO) and observational Hubble data, in order to impose complete constraints on the model pa- 
^ . rameters. We find that the power-law exponent is /3 ~ — 6.511q'25j while the Big Rip is realized 

at t s ~ IO4.5I2 Gyr, in la confidence level. Providing late-time asymptotic expressions, we 
find that the dark-energy equation-of-state parameter at the Big Rip remains finite and equal to 
wde ~ —1.153, with the dark-energy density and pressure diverging. Finally, we reconstruct the 
' phantom potential. 

(N ; 

PACS numbers: 98.80.-k,95.36. 

o 

u 

Recent cosmological observations obtained by SNIa WMAP Q, SDSS [|[ and X-ray .4] indicate that the 
observable universe experiences an accelerated expansion. Although the simplest way to explain this behavior is the 
consideration of a cosmological constant [f| , the known fine-tuning problem Q led to the dark energy paradigm. The 
dynamical nature of dark energy, at least in an effective level, can originate from a variable cosmological "constant" 
Q, or from various fields, such is a canonical scalar field (quintessence) Q, a phantom field, that is a scalar field with 
a negative sign of the kinetic term [j| [l(]| , or the combination of quintessence and phantom in a unified model named 
^v^j ' quintom Finally, an interesting attempt to probe the nature of dark energy according to some basic quantum 
gravitational principles is the holographic dark energy paradigm [l2| (although the recent developments in Horava 
gravity could offer a dark energy candidate with perhaps better quantum gravitational foundations [IH). 

The advantage of phantom cosmology, either in its simple or in its quintom extension, is that it can describe the 
phantom state of the universe, that is when the dark energy equation-of-state parameter lies below the phantom divide 
— 1, as it might be the case according to observations [1]. Additionally, a usual consequence of phantom cosmology in 
its basic form, is the future Big Rip [l4| or similar sing ularities [l5j , and thus one needs additional non-conventional 
i—< • mechanism if he desires to avoid such a possibility [l6[ . 

On the other hand power-law cosmology, where the scale factor is a power of the cosmological time, proves to be 
• '■h . a very good phenomenological description of the universe evolution, since according to the value of the exponent it 
can describe the radiation epoch, the dark matter epoch, and the accelerating, dark energy epoch [l7T - fl9l ]. Although 
it is tightly constrained by nucleosynthesis [20L 1 211 . considering late universe, it is found to be consistent with age 
of high-redshift objects such as globular clusters |22j] . with the SNIa data [H, H3|, and with X-ray gas mass fraction 
measurements of galaxy clusters I25I 26 1. Furthermore, in the context of the power-law model, one can describe 
the gravitational lensing statistics 27 1, the angular size-redshift data of compact radio sources [28[, and the SNIa 
magnitude-redshift relation [2^, [27} ■ 

In this work we desire to impose observational constraints on phantom power-law cosmology, that is on the scenario 
of a phantom scalar field along with the matter fluid in which the scale factor is a power law. In particular, we 



use cosmological observations from Cosmic Microwave Background (CMB), Baryon Acoustic Oscillations (BAO) and 
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observational Hubble data (Hq), in order to impose complete constraints on the model parameters, focusing on the 
power-law exponent and on the Big Rip time. 

This paper is organized as follows. In section |n] we construct the scenario of phantom power-law cosmology. In 
section HTT1 we use observational data in order to impose constraints on the model parameters, and in section HVl we 
discuss the physical implications of the obtained results. Finally, section IVl is devoted to the conclusions. 



II. PHANTOM COSMOLOGY WITH POWER-LAW EXPANSION 



In this section we present phantom cosmology under power-law expansion. Throughout the work we consider the 
homogenous and isotropic Fricdmann-Robertson- Walker (FRW) background geometry with metric 



ds 2 



dt 2 - a 2 (t) 



dr 2 

^ + r 2 dn 2 

- kr 2 2 



(1) 



where t is the cosmic time, r is the spatial radius coordinate, D,2 is the 2-dimensional unit sphere volume, and k 
characterizes the curvature of 3-dimensional space of which k = —1, 0, 1 corresponds to open, flat and closed universe 
respectively. Finally, as usual, a(t) is the scale factor. 

The action of a universe constituted of a phantom field 0, minimally coupled to gravity, reads [lol ] : 



S = 



R 



16ttG 



,<f> + v{4>) + l, 



(2) 



where V{4>) is the phantom field potential, R the Ricci scalar and G the gravitational constant. The term L m accounts 
for the total (dark plus baryonic) matter content of the universe, which is assumed to be a barotropic fluid with energy 
density p m and pressure p m , and equation-of-state parameter w rn = p m /p m . Finally, since we focus on small redshifts 
the radiation sector is neglected, although it could be straightforwardly included. 
The Friedmann equations, in units where the speed of light is 1, write: 



H 2 
H 



8irG 



(pm + P*) ~ ^ 

-AirG\p m + p m + p<j>+p^ + 



(3) 
(4) 



where a dot denotes the derivative with respect to t and H = a/a is the Hubble parameter. In these expressions, 
and Prf, are respectively the energy density and pressure of the phantom field, which are given by: 



1 



1 



V{<t>) 
V(<f>). 



The evolution equation of the phantom field, describing its energy conservation as the universe expands, is 

P0 + ZH(p^ +p<p) = 0, 

or written equivalently in field terms: 

dV 



■ + ZH<t> 



= 0. 



(5) 
(6) 

(7) 
(8) 



Note that as we mentioned in the Introduction, in phantom cosmology the dark energy sector is attributed to the 
phantom field, that is pue = P<t> and pde = P4>, and thus its equation-of-state parameter is given by 



w DE 



Pde 
Pde 



Finally, the equations close by considering the evolution of the matter density: 

p rn + 3i? (1 + w m )p m = 0, 



(9) 



(10) 
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with straightforward solution 



PmO 



(11) 



where n = 3(1 + w m ) and p m o > is the value at present time to- 

Lastly, note that we can extract two helpful relations, namely from we obtain 



P4> 



3 V _8^ 



8ttG 



while from we acquire 



1 / • k \ n 



(12) 



(13) 



Let us now incorporate the power-law behavior of the scale factor. In the case of quintessence cosmology, the 
power-law ansatz takes the usual form 



a(t) = a Q 



(14) 



with ao the value of the scale factor at present time to- However, in the case of phantom scenario, the power-law 
ansatz must be slightly modified, in order to acquire self-consistency. In particular one rescales time as t — > t s — t, 
with t s a sufficiently positive reference time, and thus the scale factor becomes [l?], HH : 



a(t) = a 



while the Hubble parameter and its time-derivative read: 



t s -t 

ts — to 



H(t)^ a ^- = -J^ 
a(t) t s — I 



H 



(t s - tf 



(15) 

(16) 
(17) 



Therefore, for /3 < we have an accelerating (d(t) > 0) and expanding (d(t) > 0) universe, which possesses additionally 
a positive H(t) that is it exhibits super-acceleration [301 ] . That is, in phantom power-law cosmology, expansion is 
always accompanied by acceleration. Furthermore, with /3 < 0, at t — t s the scale factor and the Hubble parameter 
diverge, that is the universe results to a Big Rip. These behaviors are common in phantom cosmology [10 . |31| and 
their realization is a self-consistency test of our work. On the other hand, note that the quintessence-ansatz (fT4")l 
cannot lead to acceleration or to Big Rip and this was the reason for the introduction of the phantom power-law 
ansatz ([15]) in [l7l.[29j. 

Having introduced the power-law ansatz that is suitable for phantom cosmology, we can easily extract the time- 
dependence of the various quantities, which re-expressed as functions of the redshift can be confronted by the obser- 
vational data. In particular, substituting (ITU . (|T2"|) . (fTB")) in ([5]) we obtain 



v(4>) 




(18) 



In the following we consider as usual the matter (dark plus baryonic) component to be dust, that is w m ~ or 
equivalently n — 3. Thus, using the ansatz (fT5j). and restoring the SI units using also Mp = Hc/8irG, we obtain 



3/3 2 - 2kc 2 (t s - to) 2 ? 



(t s -tf a 2 o (t s -tf0 



Pm0 C 2 (ts-tp) 3 ? 

2 al(t s -t)W 



Additionally, solving equation (|1 3[) for the phantom field and inserting the power-law scale factor, gives 



l2M 2 c 
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kc 2 (t s -t ) 2 ? 



(t s - tf 



a 2 (t s - t)W 



PmOC 2 {t s - to) 3 ? 



dt. 



(19) 



(20) 
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Finally, the time-dependence of the phantom energy density and pressure can be extracted from ([5]) and ([6]) , using 
(HU) and (HOI), namely: 



Mi 



3f3 2 



2/3 



3kc 2 {t s - 1 ) 

2/3 



(t, - t) 2 a 2 (t s - t) 



PmQC 2 {t s — to) 



3/3 



%(t s - t) 



(21) 



3@ 2 - 3/3 



(t s - t) 2 



PmOC 2 (t s — tg) 

2a 3 (t s -t) 



3/3 



3/3 



(22) 



and thus we can straightforwardly extract the time evolution of the dark energy equation-of-state parameter through 
(O as WDE^t) — p<p(t) I p$(t). Note that at t — > t s , apart from the scale factor, p^ and p^ diverge too, however wde 
remains finite. This is exactly the Big Rip behavior according to the classification of singularities of [15|. 

All the aforementioned time-dependencies can be expressed in terms of the redshift z. In particular, since 1 + 2 = 
ao/a, in phantom power-law cosmology we have 



t = t s -(t s -t )(l + z) 



-1//3 



(23) 



Therefore, using this relation we can extract the z-dependence of all the relevant quantities of the scenario at hand, 
which can then straightforwardly be confronted by the data. 



III. OBSERVATIONAL CONSTRAINTS 



In the previous section we presented the cosmological scenario in which the dark energy sector is attributed to 
a phantom scalar field, and where the scale factor is a power law of the cosmic time. Thus, in the present section 
we can proceed to confrontation with observations. In particular, we use Cosmic Microwave Background (CMB), 
Baryon Acoustic Oscillations (BAO) and Observational Hubble Data (Ho), in order to impose constraints on the 
model parameters, and especially to the power-law exponent (3 and to the Big Rip time t s . Finally, we first obtain 
our results using only the CMB-WMAP7 data [32^ . and then we perform a combined fit using additionally the BAO 
[H and Ho ones || . 

We mention that in the present work we prefer not to use SNIa data as in the combined WMAP5+BAO+SNIa 
dataset |35j . This is because the combined WMAP5 dataset uses SNIa data from (3(| |37[ which do not include 
systematic error, and the cosmological parameters derived from the combined WMAP5 dataset also differ from those 
derived from other compilations of SNIa data [38] . Inclusion of the SNIa systematic error which is comparable to the 
its statistical error can significantly alter the value of the equation of state 39]. Furthermore, recent analysis shows 
that the value of the equation-of-state parameter derived from two different light-curve fitters could be different from 
the one derived from two different datasets. This could make it difficult to identify if wde is phantom, since its 
obtained values from the two fitters are different [40] . A very recent critics on SNIa data analysis has been presented 
in [4l|. Definitely the incorporation of SNIa data in constraining phantom cosmology is a subject that deserves further 
investigation. 

Similarly to the non-phantom case [42J, the exponent ft can be straightforwardly expressed as 

P=-H (t s -t ), (24) 

where, as usual, we use the subscript to denote the value of a quantity at present, and we moreover set ao to 1. 
Furthermore, we introduce the usual density parameter Vl m = 8irGp m / (3H 2 ), and we split Q m in its baryonic and 
cold dark matter part, Of, and Vic dm respectively (f2 TO = O/, + Vic dm)- Lastly, it proves convenient to introduce the 
critical density p c — 3H 2 /8irG, and thus we can use the relation p m o = Vl m oPcO- 

In a general, non-flat geometry the Big Rip time t s cannot be calculated, bringing a large uncertainty to the 
observational fitting. However, one could estimate it, performing some very plausible assumptions [9j. In particular, 
assuming a flat geometry, which is a very good approximation [391 ] , and assuming that at late times the phantom dark 
energy will dominate the universe, which is always the case in phantom models, t s can be expressed as @ 

ts t + ||i + wdeI^Ho 1 ^ - fl m0 )- 1/2 . (25) 

Here we have to mention that there is one last assumption in extracting this relation, namely that at late times the 
dark energy equation-of-state parameter wde approaches a constant value. Fortunately, this is always the case in flat 
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Parameter 



WMAP7+BAO+#o 



WMAP7 



to 
Ho 



13.78 ±0.11 Gyr [(4.33 ±0.04) x 10 17 sec] 



70.2+i;| km/s/Mpc 
0.0455 ± 0.0016 
0.227 ± 0.014 



13.71 ± 0.13 Gyr [(4.32 ± 0.04) x 10 17 sec] 
71.4 ±2.5 km/s/Mpc 
0.0445 ± 0.0028 
0.217 ±0.026 



TABLE I: Observational maximum likelihood values in la confidence level for the present time to, the present Hubble parameter 
Ho, the present baryon density parameter Slf,o an d the present cold dark matter density parameter JIcdmo, for WMAP7 as well 
as for the combined fitting WMAP7±BAO±Jf . The values are taken from 



Parameter WMAP7+BAO+ff WMAP7 

P -6.5llo.25 -6.5 ±0.4 

p m o (2.52 ± 0.26) x 10" 27 kg/m 3 (2.50 ± 0.30) x 10" 27 kg/m 3 

Pco (9.3lo;4) x 10~ 27 kg/m 3 (9.57 ± 0.67) x 10~ 27 kg/m 3 

t s 104.5±2.o G y r [(3- 30 ± °- 06 ) x 1Ql8 sec l 102 -3 ± 3.5 Gyr [(3.23 ± 0.11) x 10 18 sec] 



TABLE II: Derived maximum likelihood values in la confidence level for the power-law exponent f3, the present matter energy 
density value pmo, the present critical energy density value p c o and the Big Rip time t s , for WMAP7 as well as for the combined 
fitting WMAP7±BAO±H . 



power-law phantom cosmology examined in this work, as can be seen from ([2T|) . ([22]) for k — 0, recalling also that (3 
is always negative in an expanding universe. In this case, at late times we indeed have: 

w DE + (26) 

which lies always below the phantom divide as expected 1 . In addition, one can straightway extract [H(t) 2 ) through 
OH) as 

[H{t?]=H*(±^j . (28) 

Finally, as we have mentioned, the time- functions can be expressed as redshift-functions using (I23[) . 

Having all the required information, we proceed to the data fitting. For the case of the WMAP7 data alone we use 
the maximum likelihood parameter values for Hq, to, JIcdmo and flbo [39|, focusing on the flat geometry. Additionally, 
we perform a combined observational fitting, using WMAP7 data, along with Baryon Acoustic Oscillations (BAO) in 
the distribution of galaxies, and Observational Hubble Data (Hq). The details and the techniques of the construction 
are presented in the Appendix. 



IV. RESULTS AND DISCUSSIONS 



In the previous section we presented the method that allows for the confrontation of power-law phantom cosmology 
with the data. In the present section we perform such an observational fitting, presenting our results, and discussing 
their physical implications. 

First of all, in Table HI we show for completeness the maximum likelihood values for the present time to, the present 
Hubble parameter H , the present baryon density parameter fl^o and the present cold dark matter density parameter 
f^CDMo, that was used in our fitting (39|, in WMAP7 as well as in the combined fitting. In the same Table we 



1 Note that if instead of wde we consider the effective iD e ff, that is including the weighted contribution of matter, then we have 



«W -> -1 + ^ (27) 



at t — > t s [43||, for any curvature value. 
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also provide the la bounds of every parameter. In Table |TT] we present the maximum likelihood values and the la 
bounds for the derived parameters, namely the power- law exponent /3, the present matter energy density value p m o, 
the present critical energy density value p c o and the Big Rip time t s . As we observe, j3 is negative, as expected in 
consistent phantom cosmology. We mention here that the phantom power-law ansatz (|15p is technically different from 
the quintessence one (|14[) . and thus one cannot straightforwardly compare the exponent values of the two cases (for 
example a similar wde is produced by significantly different exponents in the two scenarios 23]). Now, note that 
the Big Rip time is one order of magnitude larger than the present age of the universe, which shows that such an 
outcome is unavoidable in phantom cosmology, unless one include additional mechanisms for the exit from phantom 
phase [16j, an approach that was not taken into account in this work. 

Let us discuss in more detail the values and the evolution of some quantities of interest. For the combined data 
WMAP7+BAO+H , the potential ^ is fitted as 

V (jt) « /0 6 n 47 1 X n1 1 f 7 .„ - 2.51 x 10- 371 (3.30 x 10 18 - i) 1954 , (29) 



(3.30 x 10 18 - t) 2 



while WMAP7 data alone give 



V(t) * (3.2 fx 10^-^ ~ L " X ltr368 ( 3 ' 23 X 1018 ^ *) 19 ' 39 - ( 3 °) 



Note that the second terms in these expressions, although very small at early times, they become significant at late 
times, that is close to the Big Rip. In particular, the inflection happens at 22.4^20 Gyr (WMAP7+BAO+i?o) and 
22.0+ 3 j Gyr (WMAP7), after which we obtain a rapid increase. 

Now, concerning the scalar field evolution (f>(t), at late times (t —> t s ) the p.,„o-term in (f2"0")) can be neglected. Thus, 
(l20l) reduces to 




2M| C f3 



dt, (31) 



h (t s - tf 

which can be fitted using combibed WMAP7+BAO+£To giving 

<j>{t) « -2.64 x 10 13 In (3.30 x 10 18 - t), (32) 
while for WMAP7 dataset alone we obtain 

<f>(t) fa -2.63 x 10 13 In (3.23 x 10 18 - i). (33) 

As expected, both the phantom field and its kinetic energy (— (p 2 /2) diverge at the Big Rip. 

Having fitted the phantom potential V(t) and the phantom field itself <f>(t), it is now straightforward to obtain 
the potential as a function of the phantom field, namely V(4>). In particular, (|3"2"j) and (|3"3")l can be easily inverted, 
giving t(4>), and thus substitution to (|29|) and ([30)) respectively provides V{4>). Doing so, for the combined data 
WMAP7+BAO+H the potential is fitted as 

V{4>) « 6.47 X l 27 e 0.75xl0- 1 ^ _ 2 51 x 10 -371 ^7.4x10^^^ (34) 

while for WMAP7 dataset alone we obtain 

1/(0) « 6.37 X lO 2 7 e O-76xlO- 13 _ L9g x 10 -368 e -7.4x 10^^ (35) 

In order to provide a more transparent picture, in Fig. Q] we present the corresponding plot for V(<fi), for both the 
WMAP7+BAO+i? as well as the WMAP7 case. 

Let us now consider the equation-of-state parameter for the phantom field, that is for the dark energy sector. As we 
mentioned in the end of section \H\ it is given by WDE{t) — p<j>{t)/ f><j>{t), with p,p(t) and p<j>(t) given by relations (|22|) 
and (|2"Tj) respectively. Finally, one can extract the redshift dependence using (|2"3")l . One can therefore use WMAP7 
and WMAP7+BAO+£fo observational data in order to fit the evolution of wde(z) at late times, that is for t — > t s , 
or equivalently for z —> — 1. For the WMAP7+BAO+7Jo combined dataset we find 

1 6.068 

w DE (z) « 2" 3.670-(l + z)3-307 » ( 36 ) 
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V{(f>) (J/m 3 ) 




FIG. 1: The phantom potential obtained from observational data fitting of WMAP7 and WMAP7+BAO+Ho- 
while for the WMAP7 dataset alone we have 

1 6.328 

wdbM « 2- 3 , 824 _ (1 + z)3 .30 9 - ( 37 ) 

As we observe, at t —> t s , wde becomes -1.153 for the combined dataset and -1.155 for the WMAP7 dataset alone. 
However, as we have already discussed in the end of section [TTl at t — > t s despite the finiteness of wde, the phantom 
dark energy density and pressure become infinite. These behaviors are the definition of a Big Rip 15], and this acts 
as a self-consistency test of our model. 

V. CONCLUSIONS 

In this work we investigated phantom cosmology in which the scale factor is a power law. After constructing 
the scenario, we used observational data in order to impose constraints on the model parameters, focusing on the 
power-law exponent /? and on the Big Rip time t s . 

Using the WMAP7 dataset alone, we found that the power-law exponent is j3 « —6.5 ± 0.4 while the Big Rip is 
realized at t s rs 102.3 ± 3.5 Gyr, in la confidence level. Additionally, the dark-energy equation-of-state parameter 
wde lies always below the phantom divide as expected, and at the Big Rip it remains finite and equal to -1.155. 
However, both the phantom dark-energy density and pressure diverge at the Big Rip. 

Using WMAP7+BAO+£f combined observational data we found that fi w -6.51 ±0 25, while t s « lO^^o Gyr, 
in la confidence level. Moreover, wde at the Big Rip becomes -1.153. Finally, in order to present a more transparent 
picture, we provided the reconstructed phantom potential. 

In summary, we observe that phantom power-law cosmology can be compatible with observations, exhibiting ad- 
ditionally the usual phantom features, such is the future Big Rip singularity. However, it exhibits also the known 
disadvantage that the dark-energy equation-of-state parameter lies always below the phantom divide, by construc- 
tion. In order to acquire a more realistic picture, describing also the phantom divide crossing, as it might be the case 
according to observations, one should proceed to the investigation of quintom power-law cosmology, considering apart 
from the phantom a canonical scalar field, too. Such a project is left for future investigation. 

Finally, let us make a comment on the nature of the investigated scenarios. Although the classical behavior of 
phantom fields has a very rich phenomenology and can be compatible with observations, as it is known the discussion 
about the construction of quantum field theory of phantoms is still open in the literature. For instance in [44| the 
authors reveal the causality and stability problems and the possible spontaneous breakdown of the vacuum into 
phantoms and conventional particles in four dimensions. However, on the other hand, there have also been serious 
attempts in overcoming these difficulties and construct a phantom theory consistent with the basic requirements of 
quantum field theory [45j j . with the phantom fields arising as an effective description. The present analysis is just a 
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first approach on phantom power-law cosmology. Definitely, the subject of quantization of such scenarios is open and 
needs further investigation. 
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Appendix: Observational data and constraints 

In this Appendix we briefly review the main sources of observational constraints used in this work, namely WMAP7 
Cosmic Microwave Background (CMB), Baryon Acoustic Oscillations (BAO), and Observational Hubble Data (Ho). 
In our calculations we take the total likelihood L oc e~ x / 2 to be the product of the separate likelihoods of BAO, 
CMB and H . Thus, the total x 2 is 

X 2 (Ps) = Xcmb + Xbao + Xh - (A- 1 ) 



a. CMB constraints 

We use the CMB data to impose constraints on the parameter space, following the recipe described in (35[. The 
"CMB shift parameters" [46[ are defined as: 



R — \/Q m0 H r (z») , l a = nr (z*) /r s (z») . (A. 2) 

R can be physically interpreted as a scaled distance to recombination, and l a can be interpreted as the angular scale 
of the sound horizon at recombination. r(z) is the comoving distance to redshift z defined as 

r(z) = f ^(fe, (A.3) 
Jo H(z) 

while r s (z„) is the comoving sound horizon at decoupling (redshift z„), given by 

,00 1 

r s (z*) = / dz. (A. 4) 

Jz. H(z)y/3(l + R b /(l + z)) 

The quantity Rb is the ratio of the energy density of photons to baryons, and its value can be calculated as Rb = 
31500£lboh 2 (TcMB/2-7K)~ 4 , (ilbo being the present day density parameter for baryons) using Tomb = 2.725 [35| . 
The redshift at decoupling z* (flbo, fi m o, h) can be calculated from the following fitting formula |47j |: 



z* = 1048 

with gi and gi given by: 



1 + 0.00124 (n b0 h 2 ) 0738 



l + g 1 (n m0 h 2 y 2 , (A.5) 



_ 0.0783 (n b oh 2 )~ - 238 

91 ~ 1 + 39.5 {n b oh 2 ) O W3 

0.560 

92 ~ l + 21.1(r! b o^ 2 ) 1 - 81 ' 
Finally, the % 2 contribution of the CMB reads 

Xcmb = v cMBCinvV C MB- (A. 6) 
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Here Vcmb = P — Pdata, where P is the vector (l a ,R,z*) and the vector Pdata is formed from the WMAP 5-year 
maximum likelihood values of these quantities [35| . The inverse covariance matrix C; nv is also provided in 35] . 

b. Baryon Acoustic Oscillations constraints 

In this case the measured quantity is the ratio d z = r s (zd) /Dy {z), where Dy (z) is the so called "volume distance", 
defined in terms of the angular diameter distance Da = r (z) / (1 + z) as 



D v (z) 



(l + z) A D\(z)z 
H(z) 



1/3 



(A.7) 



and Zd is the redshift of the baryon drag epoch, which can be calculated from the fitting formula |48j [: 



where b\ and 62 are given by 



Zd 



1291 (n m0 h^y- 251 

1 + (O A/0 ^) - 828 



1 + 61 (n w h 



2 \b 2 



(A. 



61 = 0.313 (n m0 h 2 y 



1 + 0.607 (fl m0 h 



2 \ 0.674 



6 2 = 0.238 (flmoh 



2 \ 0.223 



We use the two measurements of d z at redshifts z = 0.2 and z — 0.35 |33j| . We calculate the % 2 contribution of the 
BAO measurements as: 



^BAO C'inv VbAO ■ 



Xbao 



(A.9) 



Here the vector V B ao = P Pdata, with P = (do. 2, do. 35), and Pdata = (0.1905,0.1097), the two measured BAO 
data points (33j . The inverse covariance matrix is provided in (33j . 

c. Observational Hubble Data constraints 

The observational Hubble data are based on differential ages of the galaxies In Jimenez et al. obtained 
an independent estimate for the Hubble parameter using the method developed in (4t|, and used it to constrain the 
equation of state of dark energy. The Hubble parameter, depending on the differential ages as a function of the 
redshift z, can be written as 



H(z) 



1 dz 
+ z dt 



(A.10) 



Therefore, once dz/dt is known, H(z) is directly obtained [511 ] . By using the differential ages of passively-evolving 
galaxies from the Gemini Deep Deep Survey (GDDS) [13] and archival data (5^, Simon et al. obtained H{z) in the 
range of < z < 1.8 [5l|. We use the twelve observational Hubble data from |34j listed in Table Hill 



z 





0.1 


0.17 0.27 


0.4 


0.48 


0.88 


0.9 


1.30 1.43 


1.53 


1.75 


H(z) (km s^MpcT 1 ) 


74.2 


69 


83 77 


95 


97 


90 


117 


168 177 


140 


202 


la uncertainty 


±3.6 


±12 


±8 ±14 


±17 


±60 


±40 


±23 


±17 ±18 


±14 


±40 



TABLE III: The observational H(z) data [34 



The best-fit values of the model parameters from observational Hubble data 51( are determined by minimizing 



i=l 



[H t h{p s ;z l ) - H obs (zj)} 
a 2 { Zl ) 



(A.11) 
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where p s denotes the parameters contained in the model, H t h is the predicted value for the Hubble parameter, H t, s 
is the observed value, cr(zi) is the standard deviation measurement uncertainty, and the summation runs over the 12 
observational Hubble data points at redshifts z$. 
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